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ABSTRACT 

We have analyzed Chandra High Resolution Camera observations of the starburst 
galaxy M82, concentrating on the most luminous x-ray source. We find a position 
for the source of R.A. = 09'^55™50".2,decl. = +69°40'46".7 (J2000) with a la radial 
error of 0.7". The accurate x-ray position shows that the luminous source is not at 
the dynamical centre of M82 nor coincident with any suggested radio AGN candidate. 
The source is highly variable between observations, which suggests that the source is 
a compact object and not a supernova or remnant. There is no significant short term 
variability within the observations. Dynamical friction and the off-center position place 
an upper bound of 10^ — 10^ Mq on the mass of the object, depending on its age. The x- 
ray luminosity suggests a compact object mass of at least 500 Mq. Thus, the luminous 
source in M82 may represent a new class of compact object with a mass intermediate 
between those of stellar mass black hole candidates and supermassive black holes. 

Key words: black hole physics - galaxies: individual: M82 - galaxies: starburst - 
galaxies: stellar content - X-rays: galaxies 



1 INTRODUCTION 

One of the most enigmatic results to emerge from X-ray 
population studies of spiral and other luminous star forming 



galaxies is the discovery or unresolved X-ray sources which 
appear ' Lo Imvti luiiiliiUHlLlys facLuis of 10 lo lOO's Lliues 



bright X-ray source that dominates the central region of the 
nearby starburst galaxy M82. Previous Einstein, ROSAT, 
and ASCA observations have shown th at this source is vari- 
able and is close to the centre of M82 (IWatson, Stanger, & 



Griffiths 1984; CoUura et al. 1994 



Ptakfc Griffiths 1990. It 



the Ediilmgton luminosity for a neutron star (e.g. Roberts 



has b een interpreted as a low luminosity A GN ( Tsuru et al 



& War[vnik mw, zu^uy, Gyui'jJtiiiLuiJuuiuy, & wtii'd 1999; 

Wang, 



1997), a highly x-ray luminous supernova (Stevens, Strick 



mmlei, & Plelsch 1999; ColbeiL fc MuslioLzky 1999; 
Fabbiano, Schweizer, & Mackie 1997; Marston et al. 1995; 
for a review of early results see Fabbiano 1989). The origin 
of such sources is controversial. Some are located near the 
dynamical centre of the host galaxy, and hence may be low 
luminosity AGN. However, many are well outside the central 
regions of the galaxies and require an alternative explana- 
tion. Some of these highly luminous x-ray sources may be 
very luminous super nova remnants exploding into a dense 
interstellar medium ( Fabian fc Terlevich 1996| ) , or they may 
be accretion powered binary sources, in which case they are 
excellen t black hole candidate s with masses near or above 
10 Mq (Makishima et al. 2000). Deciding between these var- 



land, & Wills 199!: ), and an accreting black h ole with a mass 



ious alternatives has been complicated by the limited spatial 
resolution of pre-Chandra X-ray missions. 

One of the most extreme and controversial examples of 
a highly luminous x-ray source in a nearby galaxy is the 



in excess of 46OM0 ( [Ptak fc Griffiths 1999) ). In this paper 
we discuss early Chandra observations of M82 made using 
the High Resolution Camera (HRC; Murray et al. 1997). 
The central x-ray 'source' in M82 is resolved into several 
sources in the HRC observations. We present an analysis of 
the brightest Chandra source. Our results suggest that this 
source may be a black hole with a mass intermediate be- 
tween stellar-mass Galactic x-ray binaries and supermassive 
black holes. We describe the observations in § 2, our analysis 
in § 3, and conclude in § 4. 



2 OBSERVATIONS 

M82 was observed with the Chandra X-Ray Observatory 
(CXO; Weisskopf 1988) using the High Resolution Camera 
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Figure 1. The central region of M82 from the October 1999 observation. The contours were calculated using a counts map with 0.53" 
pixels smoothed with a Gaussian with FWHM = 1.06". The contours indicate 2.5, 5, 10, 40, and 160 counts per pixel in the smoothed 
map. The position of X41.4+60 is marked with an 'X'. Crosses indicate positions of other sources. 



(HRC; Murray et al. 1997) and the High-Resolution Mir- 
ror Assembly (HRMA; van Speybroeck et al. 1997) on 1999 
Oct 28 04:24 UT to 14:48 UT for an exposure of 36 ks and 
on 2000 Jan 20 14:51 UT to 20:25 UT for an exposure of 
18 ks. The HRC is a microchannel plate imager having very 
good spatial and time resolution, but essentially no energy 
resolution. Each photon detected by the HRC is time and 
position tagged, making possible timing studies of individ- 
ual sources in crowded fields. The HRC contains a wiring 
error, discovered after launch (Murray et al. 2000), which 
induces a 3-4 ms error in the event time tags for this obser- 
vation. As we restrict ourselves to frequencies below 1 Hz, 
this error has no effect on the analysis presented below. The 
HRC position tags have a precision of 0.132", referred to as 
'one pixel'. This resolution oversamples the Chandra point 
spread function (PSF) which has a half-power diameter of 
0.76" (Jerius et al. 2000). We used a 15.6 pixel radius to 
extract source light curves. 

We applied aspect to X-ray events from the HRC and 
filtered the data using event screening techniques (Murray et 
al. 2000) to eliminate 'ghost' events produced by the HRC 
electronics. An image for each observation was generated 
from the filtered event lists, see Fig. nl. We used the stan- 
dard Chandra software routine wavdetect to search for and 
determine the position of point sources (CIAO VI. 1 Soft- 
ware Tools Manual). We found several sources in each ob- 
servation including both transients and persistent sources. 



Here, we concentrate on the brightest source found. The 
other sources, including spectroscopy from observations with 
the Chandra Advanced CCD Imaging Spectrometer (ACIS; 
Bautz et al. 1998) , will be described in a forth-coming paper 
(|Ward et al. 200l|). 



3 RESULTS 

The brightest source in both observations is at a location 
of R.A. = 09''55""50".2,decl. = -H69°40'46".7 (J2000). Fol- 
lowing the convention of naming sources in M82 via their 
offset from R.A. = 09''51™00', decl. = -|-69°54'00" (B1950), 
we refer to this source as X41. 4-1-60 in the remainder of 
the paper. For wider use, we also denote the source as 
CXOU J095550.2-I-694047. The position uncertainty is dom- 
inated by the accuracy of the aspect reconstruction which 
we take to have a la radial error of 0.7" (Aldcroft et al. 
2000). 

The source lies 9" from the kinemat ic centre of M82 
(Weliachew, Fomalont, fc Greisen 1984), 12" from the 



2.2)1171 peak (Rieke et al. 1980), 4 from the very lumi 



nous s upernova remnan t 41.954-57.5 (Kronberg & Wilkinson 
197^; IWills et al. 1997 |), and 13" from the suggested AGN 



candidate 44.01-1-59.6 ( [ Wills et al. 1997 ; seaquist, Frayer, & 
Frail 1997J; IWills et al. 1999|). The radio source 41.31-f59 



which is likely a compact supernova remnant ( Muxlow et al 
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1994; Allen & Kronberg 1998 ) , lies near the edge o f the error 



circle. The high ly variable radio source 41.5+59.7 (Kronberg 
& Sramek 1985) lies within the error circle. This radio source 
was bright in one observation in 1981 but not detected one 
year later or sub sequently, and has bee n interpreted as due 
to a supernova (Kronberg et al. 200C). However, as only 



one detection is available and unique identification is not 
possible based on the radio spectral index alone; t he source 
may b elong to a different class of radio transient (Muxlow 



et al. 1J94). If the 1981 radio event was a supernova, it is 
likely unrelated to the x-ray source as a bright x-ray source 
was detected at thi s position in 1979 with the Einstein High 
Resolution Imager (Watson, Stanger, & Griffiths 1984). The 



408 MH z radio flux at th e position of X41. 4-1-60 is below 
2 mjy (Wills et al. 19971). The 6 cm radio flux is below 



2 mJy except during 1981 (Kronberg et al. 200C). 

X41.4+60 is a highly variable x-ray source. In the first 
observation, the HRC count rate from the source is 0.07 c/s, 
while in the second it is 0.52 c/s - a factor of 7 brighter. In 
the first observation, X41. 4-1-60 accounts for roughly 40% 
of the counts within 6", i.e. comparable to the resolution 
of the ROSAT HRI, of the source position and only 8% of 
the counts within 4', i.e. comparable to the resolution of the 
ASCA SIS. In the second observation, X41. 4-1-60 accounts 
for more than 90% of the counts within 6" and 40% within 
4'. Thus, even in the brightest states of X41.4-I-60, ASCA 
spectra of the source are significantly contaminated by flux 
from other point sources and diffuse emission. The time scale 
of the variability places an upper limit on the size of the 
emitting region of 0.08 pc. 

Using various spectral models consistent with the spec- 
trum of this source extracted from a Chandra ACIS obser- 



vation of M82 (Ward et al. 2001), we estimate that 1 c/s 
in the HRC corresponds to an observed flux of 0.9 — 1.4 x 
10"^" erg cm~^ s~^ in the 0.2-10 keV band. The range in the 
conversion factor is due to uncertainty in the ACIS spectral 
flts. The source flux in the ACIS observation corresponds to 
an HRC rate of 0.03 c/s, so the true conversion factors may 
differ from these values if the source spectrum varies with 
flux. The flux in the second observation is comparable to the 
highest fluxes observed from th e central source in M82 with 
ASCA (IPtak fc Griffiths 1999|). Ta king a distance to M82 
of 3.63 Mpc ( p^'reedman et al. 1994 ), the inferred isotropic 
source luminosity from the absorbed flux in the two observa- 
tions would then be 1.0 — 1.5 x 10' 
band for the first observation and 7 
the second. Correcting for absorption would increase these 
luminosities; conversely, the true luminosity may be lower if 
the x-rays are beamed. These luminosities are near or above 
the highest values found for non-nuclear sources in a R OSAT 
sample of nearby galaxies (Roberts & Warwick 2000). 



■''° ergs"^ in the 0.2-10 keV 
11 X 10*%rgs-^ for 



For each observation, we extracted a light curve for 
X41. 4-1-60, see Fig. EI The source appears to have roughly 
constant fiux in both observations. In particular, the light 
curve for the January 2000 observation shows no evidence 
of significant change in flux level over the observation, so 
it is unlikely that the high flux represents a flare of short 
duration. In both observations, the power spectra show no 
significant short term variability with the power comparable 
to the Poisson noise limit over the frequency range 0.0005- 
1 Hz. 

We note that oscillations with a frequency near 600 s 




Time (ks) 



Figure 2. Light curve for X41.4+60 from the HRC observation 
in January 2000. The time bin size is 500 s. Statistical error bars 
are shown. 



were reported in a preprint of this manuscript. In subsequent 
analysis, these oscillations were found to be due to instru- 
mental effects related to the HRC response to bright sources 
and the algorithms used to screen 'ghost' events produced 
by the HRC electronics. 



4 DISCUSSION 

The luminosity of X41. 4-1-60 in the January 2000 observa- 
tion, given the assumptions concerning the spectral shape 
and isotropic emission noted above, corresponds to the Ed- 
dington luminosity for a 500 — OOOMq object. The strong 
variability between observations argues against this lumi- 
nosity being due to an aggregate of sources. The variability 
also suggests that the source is a compact object and not a 
supernova, although the possibility of a supernova expand- 
ing into a highly non-uniform medium cannot be excluded. 
Soft gamma repeaters (SGRs) produce sufficient flux; how- 
ever, the longest, bright, so-called "giant", outbursts from 
SGR s last only ~ 300 s over which they show substantial de- 



cay ( [Hurley et al. 1999| ). The fact that X41.4-I-60 shows no 
evidence of decay or variability over 15 ks, see Fig. H, argues 
against it being a SGR. Origin of the high luminosity and 
variability in an accreting massive compact object appears 
plausible. 

Dynamical friction will cause massive objects orbiting 
in the stellar fiel d surrounding a galactic nucleus to spiral 
into the nucleus (Tremaine, Ostriker, & Spitzer 1975). The 
life time, t, before reaching the nucleus is related to mass 
of the object, M, its distance from the nucleus, and the 
stellar velocity dispersion. From the position given above 
for X41.4+60 and adopting a velocity dispersio n for M82 
of lOOkms"^ ( Gaffney, Lester, fc Telesco 1993| ), a rough 
upper bound can be placed on the mass of X41.4-)-60, M ^ 



10^MQ{t/10^°yr) 



If the object was formed during the 
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initial formation of M82 then i ~ 10^°yr hence M <y 10^ Af©. 
Higher masses are allowed if shorter life times are assumed, 
e.g. M ^ 10® M© for t ~ lO^yr. This may be possible if 
the object was formed recently, in a process likely to be 
distinct from that for the formation of supermassive black 
holes in galactic nuclei, or if the object was ejected from the 
nucleus in an encounter with one or two equally or more 
massive black holes. If the object was formed recently and 
outside of the nucleus, it is likely to be less massive than 
the super star cluste rs found in M82 - the most massive of 
which is 2 X 10® M© jSmith fc GaUagher 200ol [O'ConneU et 
al. 19931). Rapid formation of a compact object with mass 
greater than lOO Mn^ in the collapse of a super star cluster 



1999 



Taniguchi et 



appears possible ( Portegies Zwart et al. 
al. 200(j). 

In conclusion, the accurate x-ray position determination 
for the most luminous x-ray source in M82, made possible 
by the high angular resolution of Chandra, excl udes identi 



ficat ion with suggested radio AGN candidates ( Wills et al^ 
1997|). [The strong variability of the source argues ag ainst 



the po ssibility that it could arise from a su pernova (Ptak 
& Griff iths 1999| ; [Matsumoto fc Tsuru l"999| ) and the high 



flux places a lower bound on the compact object mass of 
500 Mq if the emission is isotropic. The low radio flux at 
the x-ray position and the displacement of the source from 
the dynamical centre of M82 argue against X41. 4-1-60 being 
a supermassive black hole similar to that seen at the centre 
of the Milky Way. The most plausible explanation for the 
object is that it is an accreting black hole with a mass of 
500 — 10^ Mq . The Chandra data strengthen this sugges- 
tion made previously on the basis of ASCA data by Ptak 
& Griffiths (199 9) and also clearly establis h that the source 
is non-nuclear (Matsumoto & Tsuru 1999). In this case, it 



would represent a new class of compact object ( Colbert & 
Mushotzky 1999) with a mass intermediate between those of 



stellar mass black hole candidates and supermassive black 
holes found in the centres of galaxies. Understanding the 



formation of such aii object ( Portegies Zwart et al. 1999 
Taniguchi et al. 2000) may provide insights into th e forma- 
tion of super-massive black holes in gal actic centers ( [Quinlan 
& Shap|iro 199c|; [Cebhardt et al. 2000|). 
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